EVIDENCE FOR HYBRIDIZATION BETWEEN TWO SYMPATRIC 
VIOLET SPECIES, VIOLA GRAHAMll AND V HOOKERIANA 

(VIOLACEAE), IN CENTRAL MEXICO 



Aurea C.Cortes-Palomec 1 

Department of Environmental and Plant Biology 

Ohio University 
Athens, Ohio 45701, U.S.A. 
ac399295@ohio.edu 



Harvey E. Ballard, Jr. 

Department of Environmental and Plant Biology 

Ohio University 
Athens, Ohio 4570 1, U.S.A. 
ballardh@ohio.edu 



ABSTRACT 

Viola hooheriana and V.grahamii (Violaceae) are closely related species in section Viola , subsection Mexicanae 
occurring in mountainous habitats across much of central Mexico. In various locations, these species grow sym- 
patrically, and the occurrence of many individuals with intermediate morphological characters has suggested 
hybridization between these species. Using a combination of morphological, ecological, and molecular evidence 
we evaluated two mixed populations of V.grahamii and V. hooheriana in mountains surrounding Lake Patzcuaro 
in the state ot Michoacan to determine the potential ol distinct hybrid morphologies, the strength of phenologi- 
cal dillerences and ecological preierences among the species and putative hybrid, and to evaluate the level of 
gene How between the two species. Our results indicate that hybridization is occurring and that hybrids are 
morphologically distinct from the parental taxa. Pre-mating isolation is weak and has favored extensive hybrid- 
ization due to largely overlapping blooming times, weak ecological isolation and absence of spatial isolation. 
Hybrids are morphologically and genetically closer to Viola graha mi i, but do show unique alleles as well. 



RESUMEN 

Viola hooheriana y V. grahamii (Violaceae) son dosespecies pertenecientes a la Subseccion Mexicanae dentro de 
la Seccion Viola del genero Viola. Estasespecies crecen simpatricamente en varias regiones en el centro del Mexico 
y la presencia de individuos con caracteristicas intermedias ha sugerido la existencia de hibridizacion entre las 
especies. Usando una combinacion de datos morlologicos, ecologicos y moleculares evaluamos dos poblaciones 
mixtas de V.grahamii y V. hooheriana en las montanas cercanas al lago de Patzcuaro en el estado de Michoacan 
en Mexico para determinar la presencia de morfologias caracteristicas a los hibridos, el grado de diferencias 
tenologicas y ecologicas entre las especies y el nivel de flujo genico entre ellas. Nuestros resultados indican la 
presencia de hibridos y que estos son morfologicamente distintos a los padres. Los mecanismos de aislamiento 
entre las especies parentales no son muy fuertes y esto ha favorecido la hibridizacion extensa entre ellas debido a 
tiempos de lloracion similares, muy poco aislamiento ecologtco y la ausencia de separacion espacial entre las 
especies. Los hibridos son morfologicamente y geneticamente mas similares a Viola grahamii, pero muestran 
aleleos que son unices a ellos. 



Hybridization is a common phenomenon in nature, as evidenced by an estimated 70,000 
natural interspecific plant hybrids worldwide (Rieseberg & Ellstrand 1993; Judd et al. 
1999), and the tact that between 16 and 37% of the plant families reported in dillcrent 
floras contain at least one hybrid taxon (Ellstrand et al. 1996). In some cases, hybridiza- 
tion can lead to the creation ol hybrid species and/or introgressants between hybrid de- 
rivatives and the parental taxa (Arnold 1992; Rieseberg 1995). The genus Viola L. 
(Violaceae), which comprises about 525-600 species (Clausen 1964; Ballard et al. 1 999), is 
well known for its taxonomic problems due to hybridization and introgression, as well as 
complex patterns of variation in individual traits (Brainerd 1924; Anderson 1954; Rus- 
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sell 1955; Ballard 1994; Krahulcova et al. 1996; Gil-ad 1998; Neuffer et al. 1999; Marcussen 
Gr Borgcn 2000; Calderon de Rzedowski 2001; Harmaja 2003). 

Analysis ol herbarium specimens Irom central Mexico has indicated that hybrid- 
ization is likely occurring among species belonging to Viola subsection Mcxicanae (Bal- 
lard 1994). Viola subsection Mexicanac W. Becker (sensu Ballard et al. 1999) is a mono- 
phyletic group consisting ol eight stoloniterous or roscttc-lorming species: Viola 
hcmslcyana G. Calderon, V.grahami i Bcnth., V. oxyodontis H.E. Ballard, V. hunri I is l l.B. Nr 
K„ V. hooheriana H.B.Nr 1<„ V. bcamanii G. Calderon, V.guatcmalcnsis W Becker, V nannei 
Polak. (Ballard Nr Sytsma 2000; Ballard et al. 1999), and an additional undescribed spe- 
cies (Ballard, per. obs.). Instances ol hybridization have been identified Irom several lo- 
calities in central Mexico occurring between Viola grahamii and lour other members ol 
the group, V.guatcnialcnsis, V. hcmslcyana, V. hookcriana and V. hum ilis. Occurrences ol 
these putative hybrids appear to be restricted however to the immediate areas where both 
parents occur sympatrically. 

Viola graha mi i (known in Mexico as "hojade pasmo”or "pcnsamientodcl ccrro”) is a 
frequent to locally abundant perennial distributed across much of Mexico. Specifically, 
it occurs in mountainous regions from northern Mexico south to northern Guatemala, 
and is one of the most widely distributed members ol the subsection Mcxicanac (Pig. 1). 
t thrives in dry to mesic sandy loam under varied forest canopy and along stream banks 
at elevations of 1950-3600 meters (Ballard 1994; Calderon de Rzedowski 2001), and often 
forms large mats connected via aboveground stolons. It produces chasmogamous flow- 
ers at the beginning ol the rainy season (June) for approximately lour weeks and eleisto- 
eamous I lowers for several months al terward (fune-December) until the onset ol the drv 
season (Cortes-Palomec 2005). 

As previously mentioned, one ol the tour species with which Viola grahamii appears 
to hybridize is V. hookcriana (known in Mexico as “violeta”). Viola hookcriana is an infre- 



quent to locally common perennial generally growing in small isolated populations across 
northern and central Mexico (Fig. 1). It grows in mesic loam under mixed deciduous and 
coniferous forest canopy at elevations between 1700-2500 meters (Ballard 1994; Calderon 
de Rzedowski 2001 ). Reproduction occurs via both chasmogamous and cleistogamous I low- 
ers and follows a similar pattern to that seen in V.grahami i (Cortes-Palomec, pers. obs.). 

Viola hookcriana and V.grahami i are divergent in several characters: growth habit, 
1 lower structure, and foliage, and can be separated without difficulty in both herbarium 
specimens and living populations (Table 1) (Ballard 1994). The putative hybrids have been 
identified in areas of overlapping distribution of the two species (Fig. 1), and are recog- 
nizable as exhibiting intermediate morphologies (Ballard 1994). 

The goal ol our study was to assess the presence and extent ol hybridization between 
V. grahamii and Vi hookcriana in natural populations, and identify intrinsic pre-mating 
isolation mechanisms that may be functioning between the two species. We examined 
morphological and phonological differences between the two species, characterized the 
ecological factors governing where individual species and putative hybrids occur, and sur- 
veyed genetic patterns to identily gene I low between species in two study sires in central 
exico where V.grahamii and V. hookcriana intermingle extensively. 




MATERIALS AND METHODS 



Site Establishment and Sampling 

Two studv sites were established in the mountains to the north of Lake Patzcuaro in the 
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A Viola grahamii 
A Viola hookeriana 
® V. grahamii x hookeriana 



Study sites 
Northern Michoacan 



Fig. 1 . Distribution of Viola grahamii, V. hookeriana, and V. grahamii x hookeriana across Mexico and northern Guatemala. Distribution 
from Ballard (unpublished data). 



municipality ot Quiroga, Michoacan, Mexico (Fig. 1). The sites were established in early 
June 2000 and observations were made weekly over a period of live weeks during the 
summer rainy season when chasmogamous 1 lower production was most extensive. The 
two specific sites were chosen due to the presence ol abundant populations ot the two 
species and their putative hybrids as identified through a survey ot the area conducted 
earlier in the spring. The first site (Site A) [19°4T57" N; 101°32'27" W] was situated on the 
lower southern slopes of Mt. Zirate, above the town of Santa Fe de la Laguna. The second 
site (Site B) [19°42' N; 101°35' Wl was located 8 km to the west of Santa Fe de la Laguna 
(Km 48 on the road from Quiroga to Zacapu (MX-15)). While both sites were located in 
superficially similar environments, they differed in the level ot disturbance and abun- 
dance of violet species. Viola grahamii was abundant in both sites, but areas ot visually 
pure V hookeriana separated from V. grahamii were found only in Site A. 

In Site A, four 20 m long transects were established. Two transects were arbitrarily placed 
through the middle of a visually “pure” population of the parental species (A-T1VG in a V. 
grahamii population and A-T4VH in a V. hookeriana population). Two additional transects 
(A-T2HY and A-T3HY) were placed through areas where putative hybrid plants inter- 
mingled with the parental taxa. Ten quadrats (0.5 m-) were randomly established along 
each transect with position determined using a random number table. In Site B the same 
procedure was followed except that only three transects were positioned resulting in a to- 
tal ot 30 quadrats (B-T1VG tor V. grahamii, B-T3HY ot mostly V. hookeriana and a mixture 
ot V. grahamii and hybrids, and B-T2HY consisting mostly ot putative hybrids). 
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[ai;i i 1 . Comparison of morphological characters used to distinguish between Viola grahamii a nd V.hookeriano 
(on the basis of Ballard 1994). The characters of hybrids are not included since they are not consistent (i.e., 
hybrids show varying combinations of traits, making them difficult to characterize). 



Character 


Viola grahamii 


Viola hookeriana 


Stolon position 


above ground 


underground 


Stolon nodes 


present 


absent 


Stipule adnation 


free or adnate 


free 


Pubescence on petioles and 


present [abundant] 


absent 


abaxial surface of leaves 






Leaf shape 


ovate-oblong to elliptic 


broadly ovate to reniform 


Length of pedice s (cm) 


3-12 


5-21 


Calyx pubesence 


cilia te 


glabrous 


Corolla color 


white 


white, sometimes purplish 


Elevation 


1 950-3600 m 


1 700-2500 m 



For transects representing “pure 







taxa at each study site, three plants that 



were not connected via stolons were randomly selected in each ol the 10 quadrats. One 
leal was removed trom each ol the three plants and dried in silica gel for later DNA ex- 
traction. A second leal, the largest ol each individual, was collected and labeled to corre- 
spond to the DNA sample. This leal was pressed and used tor subsequent morphological 
analysis. The same methodology was used for plants in the putative hybrid/mixed-taxon 
transects, except that 5 individuals were sampled per quadrat. 

Phenology and Pollinator Visitation 

In each of the 70 quadrats (40 from Site A, 30 from Site B), the number ol individuals 
present per quadrat and the number ol individuals with open chasmogamous I lowers 
were recorded weekly during the live weeks of I ield observations to interpret phenology. 
The percentage ol I lowering individuals, those bearing chasmogamous I lowers at some 
point during the study, was calculated lor each visit lor the two species and lor the puta- 
tive hybrids to allow tor a comparison ol phenology among the parental taxa and puta- 
tive hybrids at both sites. Quadrats with more than one taxon were recorded separately 
lor each taxon and analyzed as separate observations. A Kolmogorov-Smirnov two-sample 
test was performed in NCSS (1 Iintzc 1999) to compare phenological patterns among taxa 
at each ol the two sites. 

Pollinator visitation was assessed during wecklv visits to Site A, which was the onlv 




site to have clearly dil lercntiatecl populations of V. grahamii and V. hooker iana. A total ol 
40 hours ol direct observations ol the 1 lowering individuals were recorded with s 

o 

emphasis made to observe instances of pollinator movement between or among differ- 
ing taxa. The hours ol pollinator observation were equally divided between two daily 
periods relating to potential insect activity based on the times when previous observa- 
tions in the area suggested the potential pollinators would be more active, between 8:00 
and 11:00 am, and between 12:00 and 3:00 pm. 

Ecological Characterization 

Todcterminc it the taxa grow in modally dillerent microhabitats, inter whether ecologi- 
cal isolation might reduce gene I low between the species and the putative hybrids, and 
understand the local distribution ol putative hybrids relative to the parents, an ecologi- 
cal characterization of microhabitat in the 70 quadrats was conducted. Soil N, P, l< con- 
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centrationsand pH were evaluated using a LaMotte™ Combination Soil Test Kit (LaMotte, 
Chestertown, MD). Values were given in exact units lor pH while values lor N, P, and K 
were determined via a colorimetric test procedure from which concentration could be 
estimated within a set range, 5-75 ppm lor N, 5-100 ppm tor R and 50-200 ppm for K. Due 
to the inability ol obtaining exact concentration values for N, R and K, these values were 
range standardized and analyzed in rank form. 

Percent soil moisture was recorded using the gravimetric method described by Hadley 
and Levin (1967), with the modification that soil mass was determined after one week of 
air drying. Light availability as a direct means of inferring canopy closure was measured 
using a Model-C spherical densiometer (Vora 1988). Quadrats with more than one taxon 
were recorded separately for each taxon and analyzed as separate observations. Data were 
standardized and an analysis of variance was performed in NCSS to test tor differences 
among V hookeriana, V.grahamii and putative hybrids. 

Morphological Analysis 

From the pressed leaves, measurements were taken of the 0°, 30°, 60°, 90°, 120°, 150° and 
180° radii, following the procedure outlined by Ballard and Wujek (1994). The center of 
the leaf was considered to be the point along the midrib opposite the widest part of the 
leaf blade. Presence or absence of leaf pubescence on the abaxial surface of the leaf, a 
characteristic considered diagnostic for separating V. hookeriana (glabrous) and V. 
grahamii (pubescent), was also recorded. Measurements were standardized and a dis- 
criminant analysis (Manly 1994) was performed, both including and excluding visually 
identified putative hybrids to determine overlap in leaf characters among the two taxa 
and putative hybrids. The analysis was performed in NCSS using an automatic variable 
selection procedure. Measurements at 30° and 120° were excluded from the analysis due 
to colinearity of the data. The canonical scores generated from these analyses were used 
in an analysis of variance (ANOVA) with a Fisher’s LSD pairwise comparison analysis 
(Zar 1996). 

Genetic Analysis 

To test the presence of genetic differences among the parental taxa and identity potential 
gene flow within the putative hybrids, biparentally inherited nuclear markers, inter- 
simple sequence repeats (ISSRs) were used. ISSR data have proven to be highly effective 
in detecting hybridization and/or gene flow among closely related species, and to test 
the hybrid speciation hypothesis (Wolfe et al. 1998a, b; Wolfe & Randle 2001; Archibald 
et al. 2004; James & Abbott 2005). Additionally these markers are useful in studies of 
differentiation among closely related species for they arc able to utilize some of the vari- 
ability present at microsatellite loci without the need to develop species-specific 
microsatellite primers, an advantage when working with species which have not been 
well-studied genetically. 

Genomic DNA was extracted from the silica gel-dried leaves using a Wizard Genomic 
DNA purification kit (Promega, Madison, WI). Five different ISSR primers (McCauley & 
Ballard 2002) were initially screened to identify primers which would produce consis- 
tent and scorable polymorphic bands among the taxa. Two primers were selected: Wolte 
#99B KCAfoGG] and HB#15 [(GTGdGCf These were used tor amplification of ISSR prod- 
ucts in replicated 25 pi reactions consisting of 1 gg pi 1 diluted genomic DNA, 19 pi auto- 
claved distilled water, 2.4 pi lOx PCR Buffer (Gibco BRL), 2 pi MgCD (50 mM, Gibco BRL), 



2 pi of dNTP mix (10 mM, Fisher), 0.5 pi of BSA (Bovine Serum Albumin, 4 gg pH, Fisher), 
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0.25 j.il primer and 0.25 jil ot Taq polymerase (5 U pl ',Gibco BRL). The polymerase chain 
reaction was performed in a Stratagene RoboCycler 96 with hot-top (Stratagcne Inc, La 
Jolla, CA) and programmed tor 2 min. at 94°C; 40 (primer Wolfe #99b) 44 (primer 1 IB 
#15) x 30 sec. at 94°C, 45 sec. at 44°C, 1 min. 30 sec. at 72°C; 20 min at 72°C. 

PCR products were electrophoresed in a 1.3% agarose gel in 0.5 x TBE buffer with 
flanking 250 bp ladders (Gibco BRL). Gels were stained with a solution of ethidium bro- 
mide in 0.5 x TBE buffer for 20 minutes and imaged under UV light. Gel images were 
analyzed with BioMax ID image analysis software (Version 2.0.3, Eastman Kodak Com- 
pany, Rochester, NY) to identify and size fragments. Fragments were scored as present (1 ) 
or absent (0), with fragments comigrating at identical rates (± 10 base pairs, the general 
limit of resolution for agarose gels) considered equivalent. Statistics concerning I ragment 
occurrence, including total fragment number, distribution across taxa, polymorphic and 
fixed fragments, and number of fragments shared among taxa per site were calculated 
using TFPGA 1.3 (Miller 1997). 



Analysis used three discrete methods to evaluate relationships and similarity among 
the taxa and putative hybrids. UPGMA cluster analysis and Principal Coordinates Analy- 
ses (PCoA) were performed both within and between the two sites using NTSYS ver. 2 . 0 2 j 
(Applied Biostatistics Inc.). An Analysis ol Molecular Variance (AMOVA) was addition- 
ally performed in a hierarchical lashion among species and sites using WIN AMOVA 1.55 
(Excof lier 1993). All analyses were performed according to the methods described in 
McCaulev and Ballard (2002). 



RESULTS 

Phenology and Pollinator Visitation 

The parental species expressed divergent but non-significant tendencies in blooming time; 
in fact, in Site A V. hookeriana bloomed first, whereas in Site B it bloomed later (Fig. 2). 
The Kolmogorov-Smirnov two-sample test showed no significant dil lerenccs among the 
blooming times of the hybrids and the parental taxa in either site ip values = 1 in some 
pairwise comparisons) (Fig. 2). No pollinators were observed visiting violets at any time 
during this study. 

O J 

Ecology 

All three taxa inhabited modally dilferent microhabitat conditions inf luenced by light 

J J o 

and certain nutrients, in all cases where there were significant differences between the 
taxa, the hybrids grew in intermediate conditions to the parental species. The ANOVA ol 
the ecological factors revealed significant differences in light (p = 0.02), P (p = 0.04) and l< 



(p = 0.03) among V. hookeriana , V.grahami i and the putative hybrids; moisture ip = 0.32), 
pi 1 (p = 0.86) and N (p= 0.62) did not dif fer signif icantly. Viola hookeriana grew in higher 
concentrations of P, lower concentrations ot K and more shaded environments than V'. 
grahamii (Table 2). 

Leaf Morphology 

Two hundred and eighty eight individuals were classified into three groups correspond- 
ing to the two parental taxa and the putative hybrid individuals (Fig. 3). The discrimi- 
nant analysis indicated three marginally distinct groups (Fig. 4), with the GLM ANOVA 
and the Fisher’s LSD Multiple Comparison analysis indicating that they were well-sup- 
ported and significant ip < 0.05) groups. 
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Fig. 2. Histogram of phenology for flowering individuals in two sites of Viola grahamii, V. hookeriana and V. grahamii x hookeriana. 
Abundance of flowering individuals is summarized by transects. P-values depicting strength of similarities or differences among taxa 
using a Kolmogorov-Smirnov two-sample test are shown in the box.VH = V. hookeriana, VG = K grahamii, HY = hybrids. Total individu- 
als sampled are as follows (number of reproductive individuals): Site A: VH(AT1)=24(12); HY(AT2)=151 (76); HY(AT3)=232(134); 
VG(AT4)=146(77). Site B VG=35(19); HY=95(55); VH -142 (33). 
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T '.hi i 2. Means and SE for selected ecological parameters among V.hookeriana (VH ), V.grahamii (VG) and hybrids 
(HY). Superscript letters indicate a significant difference for a particular parameter among the taxa. Moisture 



an< 




are expressed as percents. R N and K values are along a scale ol 1 to 7(1- very low, 7= 




Taxon Moisture 


Light 


pH 


p 


N 


K 


VH 26(±0.1) 

VG 24( ± 1.1) 

HY 25(±0.1) 


1 6(±0.4) 
29( ± 0.3)'* 
29(±0.3) 


6.40(±0. 1 4) 
6.51 (±0. 1 2) 
6.43(±0.1 0) 


5.32(±0.49) 

3.82(±0.37)' 

4.40(±0.33) 


1 .28(2=0.1 0) 
1.14(2=0.09) 
1 .1 8(2=0.08) 


6.44(±0.1 6) ' 
7.0(±0.1 3) 1 
6.81 (±0.1 2) 



Genetic Data 

A total ot 242 Viola individuals were examined genetically with ISSRs: 33 Viola 
hookcriana in Site A and Id in Site B, 65 V.grahamii individuals in Site A and 35 in Site B, 
and 48 putative hybrid individuals in Site A and 42 in Site B. In Site A, 70 ISSR I ragments 
were resolved vs. 53 in Site B. The majority ol these were identical between the two sites, 
resulting in a total ot 75 unique and scorablc l ragments lor the study sites together. The 
relative I requcncy ol individual I ragments ranged trom less than 10% ol individuals (rare 
to ml requent) to 100% (ubiquitous). Distribution ot 1 ragments across taxa (Table 3) showed 
that some were spec i I ic toone ol the parental taxa (e.g.. live fragments were species-spe- 
cilic lor V. grab a mi /, and three were specific for V. hookcriana) whereas others were shared 
among taxa. Most ol the shared fragments occurred between V.grahamii and the puta- 
tive hybrids (eg., 17 fragments in Site B). Nine fragments were specific to the putative 




UPGMA cluster analysis (Fig. 5) showed strong separation ol the two study sites. All 
individuals diverged substantially between Site A and Site B, regardless of taxon. The 
PC 'oA coni mined these discrete site dil lcrcnces (Fig. 6A). Additional PCoA analyses were 



carried out separately on the two sites (Sire A (Fig. 613) and Sire B(Fig. 6C)f In both sites, V. 
grahamii and the putative hybrids showed the most extensive genetic overlap. In Site A, 
V. hookcriana mostly segregated from the other taxa, although some individuals clus- 



tered with putative hybrid derivatives, and others with V.grahamii. In Site B (Fig. 6C) 
twodistinct groups of V. hookcriana were evident, with some putative hybrids occurring 
within each ol them. In all cases, putative hybrid individuals were widely placed and 
substantially overlapping in distribution with the parental taxa, together expressing a 
greater level ol generic diversity. 

I he AMOVA (Table 4) showed that most of the variation was within morphologi- 
cally delined taxa (69%), although significant differences were detected between sites 
and among taxa. Variation between sites accounted lor 27% of the total variance. Differ- 
ences among parental taxa and hybrids accounted for only 3.72% of the variance. The 
same pattern was observed when the two sites were analyzed separately. In Site A, 94.5% 
ol the variance was due to variation within the taxa; 95.3% in Site B. If the hybrids were 
eliminated, variation was slightly reduced within taxa (90. 32% and 89.17%); however, lev- 

C/ / 

elsol variation within the species themselves remained high. 



DISCUSSION 

Our results show conclusively that hybridization is occurring between V. crahamii and 

y / O o 

V. hookcriana. Phenctic analysis of leal shape characters in the species, more exploratory 
than exhaustive, showed separation among V. hookcriana, V. grahamii and the hybrids. 
1 1 v brief individuals formed a distinct group from the parental taxa suggesting that leaf 
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Fig. 3. Representative leaves of Viola grahamii V. grahamii [.i hookeriana, and V. hookeriana. 
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Fig. 4. Results of discriminant analysis of leaf shape from individuals of Viola grahamii (white dots), V. hookeriana (black dots) and V. 
grahamii x hookeriana (gray dots). 
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Tabu 3. Number of ISSR bands found per site per taxon (V. hookeriana.V. grahamii and hybrids), as we 
ber of bands shared among the taxa. 



as num 



Number of bands 


SITE A 


SITE B 


Viola hookeriana 


3 


1 


Viola grahamii 


5 


5 


hybrids 


3 


9 


1/ hookeriana and hybrids 

* 


g* 


2 


V. grahamii and hybrids 


1 1 


17 


All individuals 


46 


20 




r~ 

- 0.01 



~l 

0.09 






0.20 



0.30 



0.40 



Genetic distance 



Fig. 5. UPGMA dendrogram depicting relationships among taxa and transects. Similarity compared using the AMOVA-derived PhiST 
distance matrix. (VHA = Viola hookeriana in Site A, HYA = hybrids in Site A, VGA = V. grahamii in Site A; the same format for Site B). 



traits are useful in distinguishing the hybrids trorn the pure species, similar to the find- 
ings seen in other Vio/ct (Russell 1954;Jonsell et al. 2000). While these dif lerences appear 
to follow a clear pattern, our potential discovery ol introgressant individuals may actu- 
ally blur the distinction to some decree. Extensive gene How was confirmed, and lias 



proceeded to the extent that the two sites harbored substantially different genotypic com- 
binations. At both sites, numerous morphologically "pure" parental individuals of the 
two species were demonstrated to be cryptic hybrid derivatives, with a preponderance 
of these resembling Verahciniii. 1 ligh levels of genetic variability within taxa mav be the 

O O O O / / 

product ol long-continued gene I low. The results from AMOVA suggest that around 95% 
of the genetic variation is due to variability within each taxon. While this high level of 
variation could be partly artifactual, owing to rapid evolution of ISSR primer sites, it is 
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Fig. 6. A) Principal Coordinate Analysis (PCoA) of genetic data using the Jaccard similarity coefficient for individuals collected in both 
sites (A and B). B) PCoA of individuals collected in Site A. C) PCoA of individuals collected in Site B. Percentages on each axis indicate 
percentage of total variation explained by that axis. Viola grahamii (white dots), V. hookeriana (black dots) and V. grahamii x hookeriana 
(gray dots). 
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Tabu 4. Analysis of Molecular Variance (AM OVA) of 242 individuals collected from both sites (A and B), based on 
2000 permutations. Mean squares (Ml )S), variance component and percentage of variance as well as P values 
are presented. A nested design was performed. 



Source of variation 


df 


MDS 


Variance component 


Percentage 


P 


Between sites 




286.41 


2.30 


26.79 


< 0.0005 


Among taxa 


4 


1 8.38 


0.32 


3.72 


< 0.0005 


Within taxa 


236 

241 


5.98 


5.98 


69.49 


< 0.0005 





more likely that hybridization and subsequent gene How have proceeded to the point 
that many morphologically “pure" parental individuals are actually productsol extensive 
gene I low, bridging the gap between parents and hybrids and obscuring taxon boundaries. 
Additional genetic evidence suggesting that the hybridization between V. erahami i 

* Law J j 

and V hookcriuna has likely been occurring over a long period of time is the presence of 

ents unique to the hybrids. It has been suggested by various authors that later- 
generation back crosses might contain only a small number of the fragments from the 

L / L ' J Q 

parental species (O’Hanlon er al. 1999). If this is the case, the development of fragments 
that are “speed ic" to the hybrids arc best interpreted as the result of later generation re- 
combination and long-term post-hybridization genetic reorganization, as suggested by 
O’Hanlon ct al. (1999). Studies ol inrrogression in German violets conducted by Neuffcr 

J 

et al. (1999) using RAPDs describe a similai ■ phenomenon in the V. hvinianti Rchb.— V. 
reichcnhac hiunu Jord. ex Bor. complex, but the fact that no unique hybrid bands were 

complex had a recent origin. Directional back- 
crossing has also been suggested towards V. luted subsp. sutict ica ( Willd.) Nyman in hy- 
brids ol this species and V. t ricolor L (Krahulcova et al. 1999). In this study, the presence 
ol many “hybrid-specific" bands suggested a significant time component in the forma- 
tion ol hybrid derivatives as well as differentiation in some. Our results suggest that in- 
trogression is occurring between V. grahaniii and V. hookeriuim, and that hybridization 
might have taken place thousands ol years in the past or may occur only rarely, but high 
cvclsof inrcrlcrri I ity and appropriate site conditions have allowed the hybrids sufficient 
time to backcross with the parental taxa. 

Hybridization has likely been lavored by the absence ol a temporal reproductive iso- 
lation between the two species. Despite a minor, non significant, shift in blooming times 
ol the two species (between populations) they have greatly overlapping blooming times, 
which provide extended opportunities for inter-taxon crossing. Thissmall shift in bloom- 
ing may be related to differences in microsite characteristics and lcvclsof human pertur- 
bation such as wood cutting for local uses, which were more pronounced in Sire B than in 
Site A. The dillcrcnces in canopy openness between the two sires due to human distur- 
bance can also explain why V. hooheriana was not common in Site B, while V.grahdmii 
was very abundant. Open environments with more light favor V.grahamii , which grow 
best in open-canopv lorest or open sites while VI hookerianu is most abundant in more 
shaded environments, mainly under the lorest canopy. Other differences were found in 
terms ol soil requirements; V. hookcridnci grows in higher concentrations ol Band lower 
concentrations of K than does V.grahamii. Soil nutrients in the area however are highly 
heterogenous (Cortes-Palomcc 2005), and while they may be affecting the distribution 



s 



/ 



CORTES-PALOMEC AND BALLARD, HYBRIDIZATION BETWEEN TWO VIOLA SPECIES IN MEXICO 



1131 



ol V hookeriana and V. grahamii, the hybrids appear to survive in various soil condi- 
tions, accounting for their wider local distribution in each site, especially in areas with 
moderate to heavy levels of disturbance. It is commonly suggested that since hybrid in- 
dividuals have partial genotypes from both parental species they can potentially grow in 
both parental habitats as well as in their own “intermediate” niche and even lead to the 
extinction of the parental taxa if one of them is rare (O'Hanlon et al. 1999; Levin et al. 
1996). In fact, hybrids between VI tricolor and V. lutea subsp.sudeticahave been suggested 
to have a wider range ol distribution than both parental species (Krahaulcova et al. 1996). 
In addition, disturbance can lead to "hybridization of the habitat” (Anderson 1954; 
Rieseberg & Ellstrand 1993), resulting in a mosaic ol microhabitat conditions that favor 
hybrids and in extreme cases may cause the local extinction of one or both of the paren- 
tal taxa (Rieseberg & Ellstrand 1993; Rhyme & Simberloff 1996). Hybrids in this case can 
be compared to alien plants in the sense that they are new elements in the f lora that can 
threaten the original flora of the region (Neulcr et al. 1999). Perturbation was high in 
both study sites, mainly due to human activity (i.c., livestock grazing and timber har- 
vest). If we consider our sampling of the study sites as representative lor the region, then 
the hybrids may be said to be as abundant or even more abundant locally in the zone of 
sympatry than the parental taxa. The heavy genetic overlap ol hybrids with V.grahamii 
may threaten that species with local genetic swamping, but since it is widely distributed 
and occurs as pure populations across large regions ol Mexico and northern Guatemala 
that is probably not a problem. Given the large number of hybrid individuals between 
the two species at the sites, and the comparatively small number of V. hookeriana plants 
at any given site, however, the latter species may indeed become imperiled (Levin et al. 

1996). 

Sharing of similar flowering times has likely favored pollen movement between the 
two species, however pollinator activity would play a key role in successfully moving the 
pollen. During the year this study took place, several potential pollinators were seen in 
the area of study, however none of them were observed visiting Viola, a finding which 
seems to be abnormal. Similar observations on blooming individuals ol one of the spe- 
cies (V. grahamii) over the next two subsequent summers did however reveal the pres- 
ence of the golden banded skipper ( Autochton cell us Boisduval and Le Conte, Hesperiidae) 
and an undescribed species of bee of the genus Dianthidium ( Adanthidium ) in the 
Megachilidae (Griswold, T. pers. com.). The skipper does not seem to be very specific and 
visits several species showing white f lowers (Cortes-Palomec 2005), but the bee does ap- 
pear to only frequent Viola and it could be favoring gene flow between these two species. 
More specific studies of the pollinators would be needed to better understand their be- 
havior. 

In conclusion, we have shown that hybridization occurs between Viola grahamii and 
V. hookeriana, and that the hybrids are morphologically distinct and intergrading from 
the parental taxa. Genetically, despite a larger similarity to V.grahamii, the hybrids have 
unique alleles. Ecologically, hybrids have a wider distribution and environmental toler- 
ance than any of the parental taxa, and at least in the area of study they arc much more 
abundant than VI hookeriana. Similar flowering phenologies and effective pollen move- 
ment between the species has likely favored hybridization over a long period of time 



between the two species. 
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